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1. INTRODUCTION 

Active magnetic bearing (AMB) has been developed in recent few decades to replace the conventional 
mechanical bearings. The main concept of an AMB is to use an electromagnetic force to support a body with- 
out any mechanical contact. In comparison with the traditional mechanical bearings with many drawbacks [1] 
and [2], AMB exhibits many advantages such as: frictionless, lubricant-free operation, active vibration control 
and unbalance compensation ability. Hence, AMB is a promising solution for high-speed applications such as 
high-speed motor [3], flywheel energy storage systems (FESS), .etc. For the proper operation of an AMB, rotor 
positioning control is a challenging task due to the inherent instability and nonlinearity of the system as men- 
tioned in [4], [5], and [6]. In order to deal with system nonlinearity, authors of [7] propose a linearized AMB 
model and employed conventional PID control, similar approach can be found in [8]. In addition to linear control 
trend, an off-line tunning technique for centralized and decentralized linear AMB controllers are demonstrated 
in [9]. When the variation of the rotor position is small, linearizing the model in a small region around an equi- 
librium point in order to use the linear control technique is an appropriate approach to stabilize the AMB system. 
However, the performance of the positioning system may degrade significantly when the operating point is far 
from the desired equilibrium point. Hence, nonlinear control technique has been studied to further improve the 
systems performance. 

In order to stabilize rotor when facing with a vibrating base, [10] design a sliding mode scheme for an 
AMB system. The closed-loop system exhibits robustness to uncertainties and external disturbances. Based on 
Lyapunov's direct method and the singular perturbation order-reduction technique, the state feedback control is 
derived in [11]. Through a set of numerical simulations, the state feedback control shows its strength in compar- 
ison with output feedback and deadbeat controls. Inspired by loop shaping properties of Ho optimization and 
disturbances rejection ability of the disturbance observer-based controller, [12] successfully develop a hybrid 
controller whose effectiveness are verified via simulations and real-time experiments. By investigating strong 
nonlinearity of 3-poles AMB, a feedback linearization law is composed in [13]. Rotor responses also are re- 
stricted to a safe distance from stator boundary, however, the requirement of initial conditions might be difficult 
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to archive in practice. In a circumstance of the AMB driving rotor with unknown mass imbalance, [14] proposes 
an adaptive mechanism to tackle this situation. Taking voltage input saturation into account, a robust fuzzy con- 
trol that is able to stabilize the AMB rotor is developed in [15]. Apart from conventional AMB, [16] suggests a 
design and decoupling control for axial AMB by combining optimal algorithm and feed-forward control. 

The paper proposes an nonlinear control approach to stabilizing AMB rotor problem based on backstep- 
ping control. A speed observer is utilized in this research to acquire rotor speed information. Barrier Lyapunov 
candidate function is embedded into design process to drive the AMB rotor trajectory inside a defined range. 
The stability of the closed-loop system is proven. A set of numerical simulations are given to verify the control 
effectiveness. 


2. MATHEMATICAL MODEL OF AMB 
Figure 1 shows the one degree of freedom (1DOF) AMB system considered in this paper. In this figure, 
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Figure 1. An AMB system. 


Fy, uj and i; are the electromagnetic force, applied input voltage and current of the corresponding coil k, 
respectively, with k = 1,2; Fy is the disturbance force caused by rotor mass and other uncertainties. Suppose 
that the displacement x of the rotor from the nominal position zo caused by the initial voltage and current 
(uo, ig). Denote x, and x2 are the air gaps between the rotor and the left and right side stators, it yields 


ti =r =t, T= ro tT (1) 
41 =to— ti, d2— 319-441 (2) 
uj-—ug—u, Uz = ug +u (3) 


By using (1)-(3), a fundamental operation gives 


i (4) 
ta = A |—Riü Ja aj Yt + ui 
di = B Riz 4 Ilao Fr) 002 U2 
where, 
u 2(xo — x) B= 2(xo + x) Ku NA (5) 
2L.(%9 —x) + K’ 2L,(ao + 2) + K’ ~ Hs 7 


In (5), ug is the permeability of air, N is the number of turns in each coil and A, is the cross-section area 
of the electromagnet. Due to the fact that only the position of the rotor and the two current 7; and t2 of the 
corresponding coils are available for measurement, a speed observer is needed to estimate the rotor speed. In 
this research, the nonlinear observer proposed in [17] is utilized for speed estimation as follows 

F41 


B= tO kr. — (6) 
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where, 
jj = —kn — kx (7) 
: i? i2 
À= —kA4 : : 8 
(xo — x)? (ao + x)? (8) 
K 
"dm Á 
where k is a positive constant. Suppose that the initial conditions of the system is 
n(0)=0, A(0)=0 (10) 
Denote the speed estimation error e as 
€—v—ó (11) 
Then, by differentiating both side of (11) and based on (4) and (6)-(9), it gives 
¿= ù- Ô 
; 2 : 2 
F f 
=0() -0 () +2 -1-0 -x 
Lo — T zo +x m (12) 


- (1+0 «kp id v) 
km 


= —ke 


Relation (12) means that the speed estimation error e exponentially converges to zero and the convergence speed 
depends on k. By using the speed observer (6), the state space model of the AMB system can be rewritten as 


6 - € 
— n4 0A ky +i" 


: C uM n VAZ 
Karehana] o(a) 
7) 


= —kn — ka (13) 
i = A Rii w =a) (ô €) uUi 
ig = B Riz | CERE (ô €)io U2 


3. CONTROL DESIGN 


It can be observed from (13) that the model of the AMB has the form of a strict-feedback system. 
Hence, backstepping control technique is chosen to stabilize the system. The control design is implemented 


step-by-step as follows. 


Step 1: The goal of the control design in this step is to drive the rotor to a desired position y, with 


minimized tracking error. Denote z; as the position error 


A =E- Yr (14) 
Differentiating both side of (14) gives 
A -—tv—9r 
=n+0At kla y) Te Ür ea 
Consider the following barrier Lyapunov candidate function 
2 
vi = gos (15) | mr (16) 


with dı > 0. By using the proposed barrier Lyapunov function, the position error z; is kept in a region restricted 
by (—ky, ky) where ky is a positive number. This means the overshoot of the AMB in transient state can be 
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handle by ky. This property is meaningful in practice since the mechanical contact between the rotor and stator 
can be avoided. Differentiating both sides of (16), it gives 


2121 e 
ay 2 


V 


. (17) 
cO RACE kzı +e Ye EY.) e? 
i k? — 22 dy 
in order to render Vi X 0, the virtual control action A, is selected as 
1 d . 
Aye ete e oie e nec nee ane (18) 


with c4 > 0. 
Step 2: As seen in the previous step, if A = A,, then V; < 0 which results in e = 0 and z; = 0. Hence, 
the goal of this step is to guarantee that \ approaches to A,. Denote z2 as 


z2 =A-— ùv (19) 


The candidate Lyapunov function is chosen as 


1 k2 RN PE 
V = zl eS) +5234 | A 
"TOC (a 2) 2/2 OK (=. z) c9) 


with də > 0. Differentiating both sides of (20) gives 


V2 = | 21 
2 Par 2222 — € (+ z) (21) 
Taking A from (17) and substitute it into (15) give 
: 2 2 dızı 
21 = O29 — c1 21 (ký — 27) 33 He (22) 
b “1 


Differentiating both sides of (19) gives 


(23) 


jj = ô +e- Ý, 
i Oz, is On 0 
in which, 
QA 1 di(k2 + z2) 
B= =| kerk? + 8ei73 = So 24 
5. ca [oh aed Ts en 
OX» 1 
d 2 
» 7 (25) 
Then, by substituting (22) and (23) into (21), it results in 
2 
: Zl € əf 1 3 
V = } 
E cu. (z mpy x) T (z =) 
Oz >. ðM“. ð`. Ý kÝ, 
| FA | 26 
12 (ges On | Oz 5 9 0 MS 
From (26), to guarantee that Va « 0, the virtual control signal Ày is chosen as 
. ð`. ` ; ðM’ — 6n YX kÝ, 
Ay = | "n Ô — Y, d | 
ks On s Oz (0 ) ie zı k? — 22 0 0 
with co > 0. If 
A= Av (27) 
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then Vo < 0. However, the solution for (27) can not be obtained since it depends on two variables 7; and 25. In 
oder to solve the abovementioned problem and to save the energy, a switching control strategy is employed in 
which either 71 or is is used. If x < 0, coil 2 is turned off which means 72 = 0, then 


i2 


À 2 —kA4 i (28) 


Substitute (28) into (27), the desired current of coil 1 is 


ij = iid = To — XX Àq + kA (29) 


Denote a = Ma + kA, it gives 


4 = ild = £o —tVa (30) 
Similarly, if x > 0, then ¿1 = 0 and hence 
À = —kà i (31) 
(zo + x)? 
The reference current for coil 2 is 
ig = i24 = £o zvV—ao (2) 


Remark: Equations (28), (31), and the definition of a indicate that square root operations in (29) and (32) 
always hold. 

Step 3: The goal of this final step is to obtain u and uz which drives 7, and tə to their desired values 
71,4 and to 4, respectively. Case 1: x < 0. Denote z3 as 


Z3 =4— lg (33) 
Differentiate (33), it gives 
33 =i — tia 
— — ARi, 4 ee T (6 4- e) d 


+ Au 


Qi d: Qi dz Qi dis | AKi Oii J 3 
AN bey Ldi phd ar 34 
On OY, 0Y, 2( TE 


in which 


o Va + (zo )3 Ja 82 (35) 
: X = (zo - D NE (36) 
Ft — (m - a) (37) 
ae = -Va + (a ne A (38) 
aie T ary ii 
ði 1 ða 40) 
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with, 
2 
ða = Ory V da co | OA, | k (41) 
On Oz 0 0 Oz 0 
Ja OAy V? 3v 
= tk 42 
ER ($=) d2 C2 k ros (42) 
On -k 9X V. f -EYX a OAs 
= -d f k 43 
OY, o( ) (m) (=) Ox E) 
da 1 Oe VE. 20307 oh 
OY, E da (2) 0 Om 0 een 
ða 1 
= => 45 
a (45) 
da Pry. OV La (98 p2 
Oz, 022 EAE Oz, ? 0z3 Lor 7T 
k? Ory 02A 0(k2 + 27) 
— + 2do(Ay —— i 4 
Y 0 15208 Auc Oz, Oz? (k2 — 22)? Ge) 
9A, 1 22,d, (3kj + 27) 
TE 47 
a a ea- ea xi 
The third Lyapunov candidate function is chosen as 
l5 e 
= - 4 
V3 Vo + 9^3 + 2kds ( 8) 


with d3 > 0. The differentiation of (48) is 


; . 3 3 1 z eN? Ory e \? 
E Nagase dt EUR a) ^ (3 n) Cae 


(49) 
To make V3 < 0, Z3 is chose as 
23 = —C€323 — d3z3F? 4 eF (50) 
where Ax ái 
F- - = (51) 


and c3 > 0. The control signal u; which stabilizes the AMB system around the equilibrium point xo in Case 1 
can be obtained from (34) and (50) as follows 


1 . ð . óði 
ui z [Amis Fo — Y, — ca23 dzzz F? | D. Àd x | 
1 Aki ði . Olde: Oiid 
! ege Uy qu ey pe. (52) 
A | X 2(xo — x)? oY, Y, oY, 
Case 2: In this case, 
r>0, i$—0,uj—0, 2¢=%0+2V—a (53) 
Similar to Case 1, the control action us in this case is 
1 : is y. 
uo =a L Gô — Y, — caz4 — deer + a 1 
1 [oi BKi ði . Oi» d. 1 lÓ0i»4.. 
| d 22 12.d yd dy | yr (54) 
B | On 2(zo +r)? OY, OY, B | ƏY, 
with c4 > 0 and 3 
BKi j 
G= 2 d (55) 


2(xo +r)? Oz 
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Similar to the control design section, the stability analysis of the AMB control system is also divided 


into two cases. 

Case I: x < 0 and t2 = 0. The closed-loop system can be described by 
dizi l 
k2—22 | 


0zi 
2 2 
ky 24 


€ 


Orv 
Oz € 


2 
zi) 
Lll Ody 2 
Z2 = —C2 Z2 — d2z» əz 


23 = i1 — tld 


Z = Oza = cya (k? 


where, 
K(6+ €) 
2(xo — x)? 


Oa. , Ora; 


i = A|-Ri 


-an| 


iiad. 


i = T À | 
E E c E T 


r—z TY, 


n | 


The velocity observer in this case is 


with, 
i ù = —kn — k?x 
. -2 
A= —kA— Gora? 


and 


AMN? aM 
! d2 (2) i baz 


k? OX, 0? X, 


0(k2 + 22) 


221d13k7 + p 
cm 
Oa 
2/a Oz 

1 Oa 
2/a Or 

1 Oa 
2/a 0n 

1 Oa 
2 /a OY, 


] 1 
ina c (zo — x) 


ay ced) 


ca ci greek) 
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(58) 


(59) 


(60) 
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(62) 


(63) 


(64) 


(65) 


(66) 


(67) 


(68) 


(69) 


(70) 
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Oii d 1 Oa 

= =f) —— — 73 
OY, (xo z) 2/a Y, ( ) 
27 d 1 Oa 
OY, (xo x) 2/a Y, ( ) 


To show the stability of the AMB system in this case, consider the following candidate Lyapunov function 


1 kg $a docs W d 
zl t | 59 0 5253 E 7 
Vei = g log (gs) JO Oe Oe (+ d; =) (79) 


Differentiate both sides of (75) yields 


2 
" Z1 € 
Voi = — 622 — 6222 — e322 — dy t m" x) 
On eV EM 3 3 3 
d E ds | z3F | g 7 
: E x) ? (s x) E Ado im) à (76) 


Since c1, co, c3, d1, d2, d3 are positive constants, then Voi < 0 which means system (56) is stable. Proof of 
bounded output is given in [18]. 


Case 2: x > 0 and tı = 0. In this case, the stability analysis can be treated in the same manner as in 
Case 1. 


5. SIMULATION RESULTS 


The closed-loop system is numerically tested to verify the ability of the proposed control design. The 
AMB used in the numerical simulation has the following parameters: 


Parameters Value 
Nominal air-gap 0.001m 
Number of turns 400 

Coil resistance 10 
Cross-section area 0.001m? 
Rotor mass 2.6Kg 
Rotor initial position | 0.0004Kg 
Air-gap permeability | 1.256.10~° 


Assume, initially the rotor is at a distance of 0.0004m away from the zo ie. attached to coil 1. The simulation is 
carried out in two scenarios. First, a conventional direct Lyapunov is applied and then Lyapunov function with 
appended term to limit system output response is embedded into the system. 


g x10" 


Displacement (m) 


Figure 2. Displacement. 
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Vactual | 


Velocity (m/s) 


0 0.05 0.1 015 02 025 03 035 04 045 05 
t(s) 


Figure 3. Actual and estimated velocity. 
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Figure 4. Current. 
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Figure 5. Control voltage. 


It can be seen from the system responses that the rotor displacement has a maximum value of 0.001m. 
The phenomenon implies the rotor hits coil 2 during transient period , this is undesirable in practice. 
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Displacement (m) 
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Figure 6. Displacement. 
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Figure 7. Actual and estimated velocity. 
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Figure 8. Current. 
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Voltage (A) 
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Figure 9. Control voltage. 


In order to maintain the rotor position in safe range, the barrier Lyapunov candidate function is em- 
ployed instead of the conventional one. System responses indicate better tracking performance. The rotor is 
kept in the predefined zone and possible mechanical contact with stator is avoided. In addition, the closed-loop 
system dynamics also prove the ability of the speed observer. It is noted that the control current and voltage are 
in a practical area. 


6. CONCLUSION 

This paper presents the positioning control design for an AMB system. To handle the inherent insta- 
bility and nonlinearity of the AMB, back stepping control technique is employed. Besides, by introducing a 
barrier Lyapunov function in control design, the mechanical contact which may happen in transient response can 
be avoided. An nonlinear velocity is also used to reduce the cost of the system. The feasibility and effective- 
ness of the proposed control strategy are verified by simulations results. The simulations show good tracking 
performances of the system and all control inputs are in applicable values. 
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